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8UUUABY 



By means of a general theorem founded on the haeis of 
Prandtlfa theory of the lifting linei a method is derived 
for obtaining the minimum induced drag of airfoils in the 
proximity of ideal internal houndariee. She theorem is 
applied to the case of an Ideal «lng->fuBelage oomhlnatlon 
consisting of a lifting line intersecting an infinitely 
long circular cylindrical fuselage to determine the effect 
of iring height on the minimum Induced drag. She case of 
ideal oombinations with constant circulation is also con- 
sidered in detail, as it has been treated erroneously in a 
previous analysis. The analysis presented here incidental- 
ly reveals some errors in previous work on aerodynamic 
theory* 

INSBOOUOSION 



The approximations of the Prandtl theory of the finite 
airfoil (reference 1) permit the general solution of the 
problem of minimum induced drag of Isolated airfoils and 
airfoil systems (references 2, 3, and 4). This problem has 
also been solved by Lennertz (reference 5) for a particular 
case of wing-fuselage interference, in which the ideal fuse- 
lage consists of an infinitely long circular cylinder with 
the airfoil in the midwlng position. These solutions sug- 
gested a theorem that would solve the problem of mlnlmam 
Induced drag for the most general case of wing- fuselage in- 
terference, in which any number of wings of any front ele- 
vation and any number of ideal fuselages (infinitely long 
cylinders) of any cross section are admitted. 

This note presents this general solution as well as 
an important application, the determination of the minimum 
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induced drag of wing— fuselage monoplane combinationB with 
ideal circular fuselages and varying wing height. In or- 
der to prove the theorem, it was found necessary first to 
repeat in a new form certain portions of the basic aerody- 
namic theory "because of an error discovered in the work of 
Trefftz (reference 3). This error does not affect Treffte's 
results hut its correction is important in the present anal- 
ysis. 

Accordingly, the first portion of this note deals with 
the derivation of analytic ezcpressions for the lift and the 
Induced drag of the finite airfoil and includes an explana- 
tion of the Trefftz error as well as that of a certain par- 
adoxical statement hy Frandtl on the application of the mo- 
mentum theorem to the flow ahout the finite airfoil." TKe 
rest of the paper contains the estahlishment of the general 
solution of the minimum drag problem and its application to 
high-wing and low-wing corn'oinat ions . including the deter- 
mination of load distributions. It was found that the in- 
terference effect for combinations with constant circula- 
tion has been treated erroneously by LennertB (reference 5). 
The corrected analysis is presented here in an appendix. 

The writer is very grateful to Dr. E. Priedrichs, Pro- 
fessor of Applied Mathematics at New lorls: University, for 
his guidance and assistance in the preparation of this note* 



LIFT AND INSUCBD DBA& 07 TBS FINITE AlfiFOIL AND 
WING-FUSELAaS COMBINATION 



For simplicity, the analytic expressions for thie~llft 
and the induced drag will be derived first for the single 
airfoil; the results will then be extended to include the 
presence^ of an ideal fuselage. 



Analytic Sxpresslon for Lift Force 

In the first approximation of Frandtl* s airfoil theory, 
the airfoil is regarded as a lifting line; that is, a lin- 
ear succession of elements of small chord, each possessing 
a certain profile and angle of attack. Weak loading is as- 
sumed and the vortex sheet produced by the motion of the 
airfoil is regarded as a semi- inf inlt e plane strip with 
straight vortex lines parallel to the direction of motion. 
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In this analysis, the lifting line is taken as lying at 
rest along the x-axla in an infinite "body of fluid that 
has a velocity of magnitude V parallel to the z-axis* at 
an infinite distance hefore the airfoil* Both the air- 
foil and the flov are assumed symmetrical with respect to 
the yz-plane. The velocity field of the fl'uid is repre- 
sented hy a vector of components, $ « ^ , V, vhere 

Q ia the velocity potential arising from the presence of 
the airfoil and its attendant vortex sheet, $x 

dS/dx, $y is d^/dy. and 3$/d z. tChe assump- 

tion of veak loading is equivalent to the inequality: 



o X 



ay 

31 < < V ^ 



(1) 



Also • at an infinite dletance hefore the airfoil. 



*x *y 



(2) 



In the application of the momentim theorem to this 
flow, the airfoil is considered enclosed in a very large 
rectangular hox of fluid with center at 0 and with faces 
A, B, C. D, S, and IP, as shown in figure 1. Ihe 
total upward force acting on the enclosed fluid is 



= - L 



A- 3 



p dz dx 



(3) 



where L is the lift on the airfoil, p is the static 
pressure of the fluid, and the auhscript A-S indicates 
that the integral extends over A with positive sign and 
over B with negative sign. The pressure is determined 
from the Bernoulli equation 



p + 



OP 



p = H - i. p - p V $^ - P ($3^^ + $^ + 



(4) 
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vhere p Is density and E is total pressure. This rela- 
tion is valid everyrrhere in the simply connected region 
outside the vortex sheet. 

Under the inequality (l)t the last term of equation 
(4) can he neglected, so that 

P = H - i p - p 7 (5) 

If this expression Is inserted in the integral of equation 
(3), It will reduce to 




p dz dx = - p 7 / / §2 dz dx 

= - p 7 



(J^ $ dx $ dx) (6) 

where i, z , s, and 4 are the edges of the hox shown in 
figure 1. If the faces 0 and D are removed to z = - oo 
and z = + CO, respectively, the velocity potential $ as- 
sumes the same constant valxie on s as on 4; and vhen the 
other four faces are removed to x = ± », y s= ± oo, 
sp-ect ively , 

J^J* p dz dx = - p 7 * dx = - p "^(^ * 

• / Odx- / $dx+ / * dx) = p 7 / ^ dx (7) 

where the added terms, / S dx and ® $ dx, vanish 

because dx vanishes on the edges 5 and e, and where the 
last integral sign and subscript denote integration around 
an infinitely large contour in the count erclockr Is e sense 

in the plane D at a = + 

The total upward force, Fy, is equal to the time 

rate of change of the y-component of the momentum of the 
enclosed fluid: 

Jy = d My/dt (8) 



■'A-B 

+ 
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As the flow is stationary, the y-oomponent is Just the 
amount transported per unit time through the sides of the 
hox, vith outgoing momentum taken positive: 



/ dx d. 



J J p $y + V) dx dy + y y ' p 

D-0 ^. ^ A-B 

J J P 4y dy dz (9) 

Of the six integrals on the right-hand side, the last four 
may be neglected together with the terms in in the 

first two, from the inequality (l). For the infinitely 
large hox, the integral taken over G vanishes as well, 
and 



1^ = p V J" J'^y ±x ±Y (10) 



where the integration extends over the entire XT plane, D, 
at z = + o», which is regarded as hounded by the trace, 
of the vortoic sheet. 

Tho analytic expression for L is ohtainod from equa- 
tions (3), (7), (8). and (10): 

L = -pTyy Oydxdy-PVjT^dx (11) 
D " 

This result can he transformed hy integration by parts: 

-pV r r *ydxdy = -PVjr ^dx+PTjT *dx (12) 

BO that H 

L = -PT(r$dx=PV /rdx (13) 

T L 

TThore L and H reprosont the left and the right odgos 
of the vortex shoot, r(x) is the distrihut ion. of circula- 
tion along the airfoil span, and where the relation, 
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has "been used, "being the velocity potential on the" up- 

per surface of the sheet and that on the lower sur- 

face. Equation (13) represents the Eutta- Joukowski la^ 
for the finite airfoil. 

In hia application of the momentum theorem, Trefftz 
(reference 3) obtained the result (13) but on tho "basis of 
two omissions whose effects cancol. Ho first omittod the 
contrlhution of tho prossuro forcos to tho lift in equa- 
tion (7) and thon omittod tho corresponding integral in tho 
partial integration of equation (11). This error ■ did not 
affect his results for he used only the form given hero hy 
equation (13)* In tho onsulng analysis, hoTTOvor, tho prop- 
er form of equation (11) Is of doclslvo importanco. 



Analytic Bxpression for Induced-Drag 7orce 

Tho analytic expression for the induced drag is ob- 
tained In a similar manner hy applying the momentum theo- 
rem to the z-component of forces, except that in this case 
the second-order terms must "be retained* Again, attention 
is first restricted to the finite box enclosing- tho air- 
foil. The force In the z-dlrection acting on the enclosed 
fluid is 

= - Di - p dx dy (15) 

J)-0 



where is the Induced-drag force acting on iHe aTrfoil, 

By the use of equation (4) this expression can be written 

as : - — - -- - 



= - '^i + 



^ J" f ($/+<&y^+<&/+37 $2) dx dy (IS) 



D-C 



The z-component of momentum transported through the sides 
of tho box per unit time is: 



^ = y y P C*55+V)^ dx dy + p (0^+V) $y dx dz 

+ / / p (*B+7) *x *y (17) 
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Canceling the terms In T°. this expression Is conveniently 
rewritten ast 



^=p(/T dx dy+ /T *ydxds+/T dyds' 

'^^L-B "^"ip-B 

+ p dy + p ^(/^ *y 

+ $y dx dz *x *y 

She quantity within the last parentheses Is, 

JJ^ In ~ JJJ" ^""^^ $ dx dy da = 0 (19) 

from the continuity of the fluid flow. In ^rhlch the surface 
Integral extends also over the airfoil surface where 3$/3n 
vanishes identically. In this expression dS Is an ele- 
ment of surface and n is the coordinate normal to the 
surface of Integration. Then from equations (16) and (18), 
the relation 

= aUg/dt (20) 

reduces to: 



Dj = i^i + *y - dx dy - p JJ" * 



, $y dx dz 



■•1 



*s *x *y (^^^ 

When all the sides of the box are removed to infinity, 

$ = $ = $ = 0 on 0 , 
* y » ^ (22) 

= 0 on D 
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so that the Integral extending over 0 ▼anlshes. While 
the last four Integrals are nominally of the same order 
as the first, heeause the sides A, B, Z* and T are now 
infinitely distant from the vortex sheet, the last four in- 
tegrals are actually of higher order and can he neglected. 



where the integration extends over the entire vertical 
plane at z = + oo hounded hy the trace of the vortex sheet. 



In order to treat the problem of ving- fuselage inter- 
ference, the fuselage is idealized in a manner due to 
Lennertz (reference 5). She fuselage is taken as an infi- 
nitely long cylinder, extending from zs=-oo to z = + ", 
of any cross section, with generators parallel to the e- 
axis. She airfoil is taken as a lifting line lying in the 
xy-plane. The vortex sheet is taken as the cylindrical sur- 
face lying between z = 0 and z = + a, passing through 
the lifting line, with generators (vortex lines) parallel 
to the 2-axis. Consequently, with such an ideal wing- 
fuselage combination, the contour bounding the plane D 
at 2 = + « consists of two parts: The cross section 

of the fuselage, denoted henceforth by the letter T and 
the trace T of the vortex sheet. The entire contour is 
designated C. 

The reason for this particular choice of fuselage 
shape is the following one. If the vortex sheet is re- 
flected in the plane, z = 0, the velocities of the re- 
sulting flow In this plane will be twice as large as those 
arising from the original vortex configuration. But the 
resulting flow is that induced by an infinitely long vortex 
sheet and must be exactly equivalent to the two-dimensional 
flow existing in the plane D at z = + <», Thus, the x- 
and the y- components of the fluid velocity in the plane z = 
0 have ;just one-half the values of the velocity components 
at the corresponding points (those with the same values of 
X and y) in the plane at g = + w. in particular, the 
downwash velocity at any point on the lifting line will be 
one-half the downwash velocity at the corresponding point 
on the trace of the vortex sheet in the plane D. For any 
other type of fuselage, the vortex lines will not be straight 



Hence, 




(23) 



Ideal 17 Ing- Fuselage Combinations 
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lines; so that in general the trace in the plane D has 
neither the same shape nor the same length as the lifting 
line and the downwash at the wing now becomes a complicated 
function of the downwash at the trace* In this way, the 
" two-dimensional** character of the problem is lost. 

* 

Evidently, the derivations of equations (11) and (23) 
for L and Bj, are unaffected by the presence of such 
ideal fuselages. Then 




(24) 



In general, neither of the last two integrals vanishes, so 
that the lift on such a wing-fuselage combination consists 
of two parts - a lift on the wings and a lift on the fuse- 
lage. The lift on the fuselage arises from the aerodynam- 
ic pressure distribution over the cylinder surface and is 
to be considered as induced by the presence of the wings, 
because the lift on an isolated fuselage of the type con- 
sidered here is zero. 7or fuselages of this type, these 
induced pressure forces are normal to the cylinder surface, 
that is, parallel to the xy-plane, and can only contribute 
to the lift. 7or any other type of fuselage, the resultant 
of the induced pressure forces on the fuselage will not be 
parallel to the xy-plane, in general, so that these forces 
will contribute to the drag of the combination as well. 



Lift and Induced Drag in Terms of a Oomplex Variable 

It is useful here to transform the expressions for L 
and in still another manner so as to employ the com- 

plex variable x + iy. The flow in the plane S Is two- 
dimensional and the stream function, ^, satisfying the 
eQ.uat ions. 



3r By * 9y " Bx 



(25) 



•r ^ 



;•>»■ 



i 
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can therefore he introduced. Then 

L = - p V J^J^ *y dx dy - p 

= p V / / ■»!/_ dx dy - D V « » dx 



i,:t But 0 



dx dy - p 

= - p V J^^lr dy - p 7 dx - '^r dy) 



As r is a rigid "boundary, it oust he a streamline of the 
flow in the plane D, so that 

>]r s const, on T 

and 



\lf dy = 0 



As the vortex sheet contains no sources or sinks, ^ is 
continuous in crossing T, so that 



dy = 0 

Thus the expression for the lift reduces to 
p^j L = 7 p ▼ / dx - \!r dy) 

tie 



= - p V R.P. / fCe) de (26) 



where the symhol R.P. represents the real part of the 
quantity following it (in general, complex); z now repre- 
sents the complex variable, x + iy, and f(8) is the 
flow function (complex potential) of the two-dimensional 
flow in the plane D, 

f(z) = $ (x,y) + i. ^ (x.y) (27) 
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She eaqpreesion for can also he written directly 

in terms of ^ or of f(z): 



■ iff I 



fi(z)I dx dy (28) 



The escpreaslona for L and in terms of f(z) 

are valid, of course* when no fuselage is present and may 
be regarded as the analogs of Blasius* formulas for the 
infinite airfoil. The results of this analysis as given 
in equations (26) and (28) are not original; both expres- 
sions have been obtained by Prandtl (reference 1) by some- 
what different considerations. 



Prandtl' B Paradox 



Prandtl (reference l) concludes .that the applicattan 
of the momentum theorem to the flow about the finite air- 
foil yields different results for the contribution of t-he 
pressure forces and the momentum transport to the lift*, 
depending on the order in which the faces of the box are 
removed to infinity. A rigorous analysis shows that this 
general conclusion is correct but that his precise state- 
ment is entirely inaccurate. In the notation of figure I. 
his statement isx 

"If an airfoil, situated in a medium unbounded in 
all directions, la enclosed in a control surface in 
the form of a parallelepiped, the application of the 
momentum theorem for steady flows yields the follow- 
ing results. If the bounding faces, A and S, C and 
D, are first removed to infinity, and then the faces 
B and T, the momentum theorem yields a momentum 
transport arising from the vortex sheet, which is equal 
to the lift. If the faces, A and B, E and 7, are 
first removed, and then the faces, C and D, the vor- 
tex sheet contributes nothing, but the momentum trans- 
port arising from the bound vortices yields the lift, 
rinally, if C and D, B and 7, are first removed, 
and then A and S, the momentum transport vanishest 
and the lift arises from the pressure forces. In 
other cases, both the pressure forces and the momentum 
transport are obtained.** 




4 



■ i 



12 



Technical Vote Vo. 812 



This statement Is inaccurate for two reaibons. The 
removal of the faces A and B to infinity does not elimi* 
nate the contrihut ion of the pressure forces* Shis error 
was also committed hy Treffts. as indicated ahove. The 
other reason is that the momentum transport « in any case* 
cannot he separated into contributions from the Tortex sheet 
and the bound ▼ortiees. Tor ezamplot considering the momen-' 

turn transport across X and T* / / p $ dy dSt 

d J J, J. 

application of the Biot-Savart law shows that arises 
from the vortex sheet and arises from both the bound 

▼ortiees and the vortex sheeta 



an 



Thus, attention can be restricted to the contributions 
of the pressure forces and the momentum transport. Trom 
the analysis presented, it is clear that once the contri- 
bution of the pressure forces has been transformed into 
the line Integral of equation (7), the faces A. B. C, S, 
and 7 can be removed to infinite distance from the air- 
foil in any order, for in the limit they contribute noth- 
ing to the lift expression. Then the expression for L 
becomes: 



J^J" $y dx dy - p 4x (29) 



L = - p Y 

where now the plane of integration B lies at any distance, 
s > 0, from the lifting line. Here the double integral 
represents the contribution of the momentum transport to 
the lift and the contour integral represents that of the 
pressure forces* 

Id order to find the ratio of these contributions, it 
is necessary to employ an exact expression for the velocity 
potential The vortex sheet is mathematically equivalent 

to a dipole layer, with dipole strength equal to the circu- 
lation. Hence, by the employment of a well-lcnown formula 
of potential theory, 

b w 

-b O Z2 

+ y^ + (s-s)^ ds dx (30) 
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vhere t is the semi span of the vlng, and r(z) Is the 
distrilsution of circulation along the span. This expres- 
sion for the potential has heen furnished "by E. Friedrichs 
and is much simpler than the Courier Integral expression 
derived "by von Earntan (reference 6). If equation (30) Is 
integrated with respect to z 

$ = I $1 + I *a 

where 

*i C^^.y) = ^ ^^^^ ^ Li^'^)^ + y°3"'' d^ (31) 
^h 

is the potontlal of the two-dlfflenslonal flow at z = + es, 
and 

h 

z 



$2 (x,y,z) = ^ y r<z) y C<2c-x)^ + y""]" 

[(x-x)^ + + z^n'^di (32) 

Evidently, 

$a ~'°»' s-s z — + OS 

I'irst, consider the plane H at an Infinite distance 
from the airfoil. Zhen $ s and 

jT $ dx = r(x) y C(x-x)^ + y^]"*^ dx dx 

Xhon 



/ 



y [(x-x) + y 3 dx — ». - it as C^c^+y^] — " 



$ dx = J P 



flo that 

$dx = ipv/ r(x)dxail, 
and from equation (29^ 




- P V / / $y dx dy = I I. 
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Now let B "be at any finite distance s from the 
airfoil, Then, 

- PV jr$dx=-|p V jT d3c - I p V jT $a 



oa 



00 



- 1 L - i PT / *a 4x 
00 



But for finite z. 



jT $g dx — as [^^ ~ 



CO 



BO that in this case, 

- P V 



DO 

// 



- P V / / $y dx dy = I L 



These distinct results yield the conclusion that if 
the faces A and B, B and I* are first removed to infin- 
ity and then the faces C and D, the pressure forces con- 
tribute one-quartor of the lift and the momentum trans- 
port contributos thr eo- quart or b , If, on tho othor hand, 
the facos 0 and S are first removed to infinity and then 
the faces A and B, E and I", in either order, the pres- 
sure forces contribute one-half the lift and the momentum- 
transport contributes one-half* In other cases, one of 
these two results vlll be obtained, depending only on 
whether S ' is tho last face to be romorod. 

Mathomatioally, this peculiar rosult arises from a dis 
continuity at z = + » in tho expression, « $(x,y,z) dx. 



Phyeically, it signifies that any distinction between the 
contributions of pressure forces and momentum transport to 
the lift is an strtificial one, at least when the airfoil 
is in an unbounded fluid* 
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She problem of minimum Induced, drag consists in minlm:- 
Islng the induced drag under the condition of given lift. 
Shis Is an isoperimetrio problem In the calculus of varia- 
tions* which consists in determining the analytic function, 
f(s)> which makes 



>x = I- ff lf*U)l ^ dx dy 



a minimum, with 

L = - p y tiz) dz 

given, or in short, 

6Di = 0, 6^Di > 0. with 6L = 0 (33) 

In the case of wing- fuselage interference, this problem 
contains mixed boundary conditions, which are conveniently 
expressed in terms of the stream function, \|/. The fuse- 
lage cross section is a rigid boundary of the flow in the 
complex z-plane, so that is constant on the cross sec- 

tion, 7; the trace of the vortex sheet contains no 
sources or sinks so that ^ is continuous In crossing the 
trace. Eenee, the boundary conditions aret 

ilf = 0 on I (34) 

^b ^ 35) 



where the subscripts a and b refer to values directly 
above and below the trace, respectively. 

This problem has been solved by Lennerts (reference 5) 
for the particular case of a midwing combination with cir- 
cular fuselage by employing the method of images which, in 
fact, is available for this one case only. By means of a 
generalization of Lennertz*s result, the solution of the 
general variational problem contained in equation (33) and 
the boundary conditions of equations (34) and (35) has been 
found to lie in the following theorem. 
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THBOBBM: The analytic function, f^t), which minim- 
izes the induced drag with gl7en lift and satisfies 
the honndary conditions. Is the sum of two analytic 
functions: one 1b the flow function of the downward 
potential flow about the fuselage houndary, the other 
is the flow function of the upward potential flow 
ahout the entire bounding contour, C, consisting of 
the fuselage cross section and the trace of the vor- 
tex sheet, where the two flows have equal and oppo- 
site velocities at Infinity. 

To prove that equation (33) is satisfied by the flow 
function, f(z)t advanced in the theorem, let 

f(z) = fi(z) + fa(z) ^ 

or 

* + ill; = (*i+*2) + 1 (il/j.+>lfa) 



(36)' 



where fi(z) = $i + is the flow function of the down- 
ward flow about the fuselage cross section, and f^(z) = 
$g + i>{rg is the flow function of the upward flow about 

the entire bounding contour C. -These functions satisfy 
the boundary conditions, for 

^1 f H^a = 0 on F (37) 

>ka = 0 o» T (38) 
and, from the regularity of fi(z) outside ?, 

ana =•- C^X °° 



where "6 Is the direction normal to the trace and pointing 
into the fluid region. 



ITow 

L = - p V ^ ( $ dx - >k dy) 
so that 

8L = - p V ^ (6$ dx - 8>k dy) (40) 



I I 
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Also, 

^± i F r C^x^ + dy 



so that 

fin, = P / / i^^ 6^^_ + b"^^) dy 



iff 
. - •//■ 

+ pyy^^^sx + ^Sy S>ky) dx dy (41) 



Using (rreen' s theorem: 



r 

av^ ds - P ^ 6^1^ ds 



0 



OB 



- f %^ (6^5 ^« - PJ^ % |i> (S^) ^8 (43) 
0 * 

where V is the direction normal to the hounding curves 
and pointing into the fluid region, and ds is a lino ole- 
mont of the bounding curvo, Sut 



^ &^ l^i. ds = 6^ ^ ds +^ 8M/ |ii 



ds 



On 7, 5^ vanishes hy virtue of equation (34) • Tho sec- 
ond integral is 

,H axl, 

ds = 0 



L 

from equation (89) • Also 



% — (6x|f) ds = 0 



0 

from equation (37), Honco, equation (42) roducos to: 
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5Di = - 9 j S^{r + ^2 (6^^) ds (43) 

Jor tlio noraal diroction cliosen horo for th.o infinitoly 
large contour, the Gaucliy-Riemann equations are: 

9$ 9^ 

— — — _ — 



ds dv 3v ds 

so that 

5Di 



= - p / [^8\[f d$i + \{fa d (6$) J (44) 



Integrating the second term by parts, 

= - 9 r {b'i/ d*i - 6$ d^g) (45) 

Co 

From the definitions of fi(z) and f3(z) given In the 
theorem, 

ff i(z) , -* icz. -cy, 
as z — ^ (46) 
j^fa(z) --1,-icz, $a — » cy, 

where c is a real, positive constant with tho dimonsions 
of volocity, Evidently, tho gonoral volocity vector, 
($3j., $y., proportional to c, so. that from the in- 

equality (l), 

c < < V (47) 

Prom an insertion of the limiting values of $1 and 
in equation (45) 




6Di = - p c «S (6$ dx - dy) 

= f SL (48 ) ■ 

from equation (40) « If the variations are restricted to 
those that mako 5L = 0, in accordance with oquatibh '(33) » 
tho function f(z), advanced in the theorem, makes 

6D1 = 0 (49) 
that is, makes "D^ an extremum. 
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In order to show that the resulting 13^ Is really a 
minimum, let 

Dl [tlr.v^'J = ^J' f ^x* + ^y') f^^ ^7 (50) 

where and \{/' are nov any two functions satisfying the 

l>oundary conditions given in equations (34) and (35). Then 

In this egtiation, let be the stream function of f(z), 

the flow function of the theorem, and 5>If any variation of 
>tf that satisfies the boundary conditions and makes 5L = 0, 
Then, from equations (41) and (49), 

2Di I^Tlf , 5\|fj = 6 »i = 0 (52) 



(51) 



80 that 



as I 5^lr,S\l' I is nonnegativo. Hence 0^ is a minimum 

for tho'-givon-^ ^^. 

A similar argument shows that this ^ (and consoquont- 
ly, f (z)) is unique to within a constant, ffor suppose 
any other stream function satisfying the boundary con- 

ditions also minimizes D^, so that 

Di J = Di [*.^ (54) 

Then the difference, 5>}r = ^' >lr, is an admissible varia- 
tion, and 

Di j = Di [^.^] + Di ^6ilf.6^1/j 

so that from equation (54) 

Di [sxjf.fit] = 0 

and 

8>1/ = - t = const, (55) 
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irhero tht last relation foXXov-« from the positive definite 
charactar of the form. Di CM'i^3 

finally. It is possible to derive a simple relation 
lintwoon the niinliaum induced drag and the (given) lift* 
noplacing 8* and S^U in equations (40) and (41) hy $/2 

and >l'/2, roapectlveXy, &L Is repXaood "by 1/2 and 

BD^ hy D^t Eonce, oquation (48) is replaced hy the rela« 



Tho theorem ostahllshod is significant in several ro- 
specte. First, it is quite general in that It applies t^ 
any comhination whatsoever, provided onXy that the fuselage 
Is of the type specified. Tor it Is clear from the mathe- 
matical analysis that the comhlnation could consist of any 
numher of fuselages, each of any cross section, of any nun- 
hor of wings, of any front elevation, and lying In different 
planes. 



Second, the theorem contains all tho previous solu- 
tions in tho prohlcm of nlnlmua Induced drag as special 
casQS, For when thoro is no fusolago, tho downward flow of 
the theorem reduces to a simple reotilinear flow and the 
upward flow Is Just that around the trace of the vortex 
shoot (of which several may be present), Ihis ease is the 
well-known condition of constant downwash derived hy Uunk 
(roforonce 2) and used hy him and others to find tho mini- 
mum drag of isolated airfoils and systems of airfoils (ref- 
erences 1 and 4). As previously mentioned, the solution 
obtained here was used by Lennertz to find the minimum drag 
of a particular ideal combination, Se Ealler (reference 7) 
has found the minimum drag of an airfoil in proximity to 
tho ground. This solution can be immediately established 
by means' of an obvious extension of the preceding theorem 
to include the presence of external boundaries. 



The theorem reduces the entire problem of minimum drag 
to tho detormination of tho roqulrod flow function, that 
is, a problon of conformal mapping, For a given typo of 
conbination, the detormination of the required mapping is 
gonorally a difficult task. Ono case of particular intor- 
ost can bo solved explicitly, nanoly, the high- or tho low- 
wing monoplane conbination with circular fusolago," Tho 
rest of this paper is llnltod to tho analysis of this case 
and some related oonsidoratiens. 



tion. 




(56) 
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HIGH- AND LOT-ffllSra M0N0FLAIT2 C0MBIKA7I01TS 
VITK CIBCULAB FUSZLAGS 



Only the high-wing coin"bination will "be treated in de- 
tail here, for, as will "be shown later, it is exactly equiv- 
alent in the th^eory to the low-wing coaTsination. The ideal 
combination is shown in figure 3, The fuselage is an infi- 
nitely long circular cylinder with axis parallel to OZ. 
The wings (lifting lines) lie along the I axis. The fuse- 
lage radius is taken as the unit length, and the semispan 
(distance from wing tip to plane of symmetry) is called h. 
Like all other lengths appearing here, "b is a nondimen- 

sional quantity. Such an ideal coahination is the first 

approximation to an actual monoplane combination with a 
long fuselage and wings of chord length small compared with 
both tho span length and tho fuselage radius. 

Tho bounding contour in tho plane at z = + •» con- 
sists of a circlo roprosenting the fuselage cross soctiom 
and two horizontal linear sogmonta (double lines) repre- 
senting the trace of the vortex sheet, as shown in figure 
3, Let p TT bo the .angle between tho positive Y axis and 
the- radius to the projection of the wing root on this 
plane, as shown. Then the height of the wings above the 
fuselage axis is cos § rr* 



The Conforraal Mapping 

In order to find the minimum induced drag in terms of 
the given lift, as well as the various related aerodynamic 
quantities for this combination, it is necessary to find 
the flow functions defined in the theorem of the preceding" 
section, Tho first of those, fi (z) , represents tho down- 
ward flow about tho fuselage contour (tho circlo of fig. 3). 
This function is: 

f, (z) = i e fz + i cos P TT r-^ r } (57) 

\ Z + i cos P TT^ 

Tho flow function, f a (z) , of tho upward flow about tho 
whole contour of figure 3 is quite complicated but can be 
obtained implicitly by conformal mapping. 

The analytic functiont 
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£ = i log ' ^ ^ (58) 

* 2 Bin p TT ® z - Bin p TT 

maps the exterior bf the whole contour in the z-plane on 
a region in the t-pla^is ■bounded "by straight lines, as shown 
in figure 4. The point at infinity in the z-plane is 
mapped on the origin of the ^-plane* and the points E and 
A of the ^-plane have the coordinates: 

£ = 1 log ; ^ + sln' g^ 

2 sin p TT .± "b - sin p tt 

This region of the t-plane is now laapped conformally 
on tho upper half of a conplox t-plano by tho Schwarz- 
Christoffel method, as shoirn in figure 5. The point G of 
tho £-plr.ne is mapped on the point at infinity of the t- 
plane , and the other two arbitrary points on the real axis 
of the t-plane are chosen as the point 0 laheled C and 
tho point +1 labeled O in figure 5. Tho differential 
equation is: 

dt (t® - n^) (t^ - 1) 

where d and n are the coordinates of B and F in the 
t-plane. Integration of equation (60) and evaluation of 
the constants yields as the mapping function. 



i = , / ^ [p log 2-±-i + (1 - p) log f±-i1 

2 sin PttL n-t t-lj 



2 sin p 

whore tho parametors d and n depend on ^ and p 
through tho relations. 



(61) 



a n Cn - p (n - D] 

d = (62) 

Cl + P (n - 1)3 



The origin of tho £-plano is nappod on a point on tho imag- 



m im 
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Inary axis of the t-plane, t = 1 sx * where sx satisfies 
the relation, 

n = 8x Qot - l) cot""* BxJ (64) 

Equations (62), (63), and (64) permit the evaluation of the 
aatheaatical paraoeters Si, n, and d in torms of the 

physical paranotors, h and ^« 

Sho f lov function in the t-plane corresponding to 
fa(2) is: 

5a (t) = fa C2(t)3 = -3-^ 5- (65) 

N« (i Sx) t^ + Bx 

where 

N(t) = e2^<*>«^" ^ _ 1 (66) 

The Uinlntuia Induced Drag 
The laininun Induced drag of the conhlnation is given 



hy 

""^nin ~ 2V 



where 

L = - P V H.P.jT f (z) da 

= - p V S.P. jT fi(z) dz 

- P V R.P. Jtaiz) dz 
Troni equation (57) *** 



- P V ^''^'J^ fi(2) dz = - Stt P 7 C (67) 

00 



Tho second integral in the ozprossion for L is ovaluatod 
hy transforming it into the corrosponding integral in the 
t-plano» 
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- P V R,P. / fs (z) dz = P V f Pa (*) (*) (^8) 

By oxpansion of (*) ^'^^L ^'C*) aliout t = i si , ono 

finds 

r 

P V R.P. 0 PaCt) 8» (fc) dt 

fain p Trf irw^lCisi)") 2 f ^" (iSi) 1 , 1 1 (69) 

= ^ b'CisiyJ iU'dsi)]" z Ln« (isiy J"" Bi^j 

T7hQre K(t) is defined by equation (65), Addin«r- ejjurttions 
(67) and (69), 



L = iTp7c/M(t.p) (70) 

irhere 



_1 ^ 2 f sin p TT l'/rK'"(iBi)l' 

M(l3.p) Lir« (isi) J \Ln' (isi) J 



Tta.en 



_ S r^l!JiSiil+ 1 "I - 2 (71) 
3 Lff' (isi) J Si^J 



^l-min " TrPc''/2M(l),p) (7.2) 



Sliminating c "betneon equatione (70) and (72) , 



Di„,„ = M(b.p) (73) 

rain 2tt P T 

Equation (7i5) ezpreaees the dopendonoo of tlie lalnimim in- 
duced drag of the corabinati on on the given lift in terms 
of the nondisensional lengths used in this section. The 
corresponding dimensional expression is: 

^i«in = ^% a M(D,p) (74) 

where R is the fuselage radius. 

Po»r given lift, varios directly with MCb,^), 

■^nin 

This quantity has hoon evaluated nunorically and is shown 



I 



XTACA Tochnlcal Noto Ho. 812 25 



in figure 6 ae a function of tbo wing hoight, cos § TTt 
for sovoral values of 'b. Sho values for negative wing 
hoiglats have 'boon found froa a rolation to "be dorivod lator 
in this Boctlon. 

Iho Intorforonco Sffoct * 

The offoct of tho presonco of the fuselage upon the 
niniiaun induced drag can "bo found 'by conparing that of tho 
conbinations xrlth that of an isolated lifting line of the 
sane span longth and total lift. Tith tho use of nondl~ 
□onslonal lengths* tho ninlnun induced drag of tho isolated 
ving is: 

4 , = 4-3- (75) 

^nin 2rr p T h 

Honoo» tho relative Incroaso in tho ninlaun induced drag 
of tho conhlnation as oonpared Tilth that of tho isolated 

iring is: 

1(1,, = ""^uln - °^nin ^ ^^^^^^ ^ ^ ^^^^ 



She dependence of this "interference coefficient" on the 
semi span h and the wing height cos p tt is shown in 
figure 7« 

The Low-ITing Combination 

She case pf the low-wing comhination is treated by 
considering a comhinatlon of semlspan b and of wing 
height - cos p tt = cos(l - tt* She minimissing flow 
function for this comhination represents the superposition 
of the downward flow about the fuselage cross section and 
tho upward flow about tho entiro contour In the z-plane. 
Vhen the axes are rotated through 180°, this oomhlnation Is 
transformed into the corresponding high-wing comhlnation of 
semlspan h and wing height cos p tt while the minimiz- 
ing flow function is transformed into - f(z), where 
f (z) is the minimizing flow function for the high-wing 
combtnationa Hence, all the relations previously obtained 
are equally valid for tho low-wing combination and, in 
particular, this argument yields the important results; 
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M(T3, p) = M(b, 1 - p) 
1(1), p) = 1(1), 1 - p) 



(77) 



These relations have already been used in plotting figures 
6 and 7* 

She complete eq.uivalenco of hlgh-\iring and loTr-wing 
comljinations in this theoretical first approximation is not 
rofloctod in oxporimental results (references 6 to 13), in 
which the presence of the ho\indary layer creates a funda- 
mental difference hetweon the ttro typos of combination. 
For unfillotod comhinations , tho oxporimonts shorr that the 
drag characteristics of high-wing comhinations are much 
superior to those of tho low-ning typo but that the lift 
characteristics are nearly tho same, tho high-vrlng ctimblna- 
tion boing only slightly superior. 



LOADING FROFERIIES 07 VIK<}-FUS2LAaS COUBIKAIIONS 



As indicated in connection trith equation (24), the 
lift on iring-fuselage combinations is composed of a lift 
force on the trings and a lift force on the fuselage. In 
this section, the distribution of these loads over the com- 
bination width is determined and, in particular, the effect 
of changes in the xring height is investigated, Bxcessivo 
calculations are avoided by treating in detail only tho 
case of tho oxtrorao high-wing combination, (p = 0) ; tho 
midwing case has already been treated by Lennertz (roforonce 
5), but some considorationt is made also of combinations with 
intormediato wing heights, This troatmont automatically in- 
cludes tho case of the extreme low-wing combination (p = l). 



Tho bounding contours in tho z-plano for such an ex- 
treme high-wing combination is shoirn in figure 8« From 
equation (24), 



Determination of Lift Distributions 



L 



Lji + Lff 



where 




I 
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is tho lift on tlxo fusolago and 



= - P 7 ^ * dac 



is tho lift on tlio wings. Writing ' 

yrhoro and aro the voloclty potentials of the down- 

vard and the upirard flows of tho thooroa, rospoctiTcly. 



^3P = ^13- ^2S 



T/horo 



Lj^j. = - p V 9& $1 dx 
Lgj. = - P T J ^3 

Th.0 di strilsut ion of lift across tho fuselage vridth is thoro- 
foro defined hy tho equation: 



dx ~ dx "** dx 



vihore dL^/dx is tho lift por unit longth in the x~direction 
and the subscripts a and h refer to the top and the hot- 
tom sides of the fuselage section. Tho quantities, dL^^p/dx 
and dLgp/dx may "bo rogardod, for the purposes of this boc- 
tion, as partial lift distributions arising from the sepa- 
rate flov7 functions, fi (z) and fgCz), From equation (57), 
Trith p = 0; 



80 that 



f,(z) = i c (z + i - -O-J) 



and 



*ia = - 2c /l - x^, = + 2c yi - 

^ = + P V ($,a - $,t,) 



2^ 

X 



- 4 P V c /l - X®, - 1 < X < 1 (78) 
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This oxprosslon represents an elliptic dlstri'bution of neg- 
ative lift across the fuselage width or, so to speak, the 
downward flow gives rise to a downward thrust on the fuse- 
lage. 

In order to find dLgp/dx* the distrihution of $a 
over the fuselage cross section must he determined* Shis 
distrihutioDi of potential is ohtalnod from the mapping 
process descrihed in tho procoding section, hy taking the 
limit p = 0* In this way, it is found that the function, 

^ = 1/z (79) 

maps the exterior of the contour in the z-plane on a region 
in the t-plane hounded hjr the straight lines shown in fig- 
ure 9, This rogioni is mapped on the upper half -plane shown 
in figure 10 hy tho function 

S = - - log g 

ttLs t+1 n^-lJ 

whore tho paramotor n d.oponds on tho somispan, h, 
through tho rolatioir., 

1 1 fl , l^ + 1 n 1 

h = - Li rrr ^rnj 



The imago in tho t-plano of tho point at infinity of the 
z-plano 11 08 on tho Imaginary axis at 



t = i Sj^ 
whore Sj. is defined hy tho relation. 



(80) 



n = /I + =1 

cot 8], 

Tho flow function in tho t-plano corrosponding to f^ (z) 
is 

2c s, 1 

I-aCt) = fa C2(t)3 = 777~~- -5-^^^ a (81) 

t* (isi) t + 8i 



ITrom equation (81), tho distribution of potontial along 
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the houndary (real axis) of the half t-plane is dotorminod, 
and through equations (79) and (80), the corresponding po- 
tential distrlhutlon along the hounding contour of the z- 
plane. Eence, the lift distrihution, 

can he found from this graphical method. Calculations have 
"been performed for the cases h = 2 and h = 6 and tho 
results are lllustratod In figures 11 and 12, which also 
show dLj^ji/dx (givon hy oquAtiort (78)) as well as tho to- 
tal lift distrihution, dLp/dz, Shoso curves show that the 
lift^ on tho fuselage, givon "by tho area under tho curvo for 
dL;gi/dx in each case, is negative as can ho demonstrated "by 

olemontary considerations • Xho curves of figuros 11 and 
12 show tho values of dLj^p/dx, dL2i>/d3C, and dLji/dx, 

each divided hy tho convoniont factor, 2tt P V c. Tho ac- 
tual valuos in each case, of course, dopond on the total 
lift on tho comhination, and may ho found from tho rela- 
tion, 

2rT P Y c = 2L/U(h,0) 

She lift distrihution on the wings is found lay tho 
same method oxcept that for tho wings only $q contrih- 

utos to tho lift, Tho results for h = 2 and h = 6 aro 
shown; in figuros 13 and 14. Thoso distributions do not 
differ markodly from elliptic distrihutions. 



Lift Distrihution ovor Tusolago Ifidth 

It has hoon shown that the lift on the fuselage is 
negative in oxtrorao high-wing comhinations with minimum 
induced drag. From elementary considerations, ono is also 
led to expect negative lift on the fuselage in similar com- 
hinations with constant circulation. This case has already 
heen investigated hy Lonnertz (roforence 5), who ohtained 
a positive fuselage lift. This result has heen found to he 
erroneous; tho corrected analysis of this case is presontod 
hero in the appendix. 

For the sako of complotonoss , tho loading distrihu- 
tions over tho coahinatlon width havo hoon plotted to a con- 
voniont soale for eight different casos in figuros 15, 16, 
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17, and 18 to show the dependence of the loading on the 
Tfing height, the span length, and tho distribution of circu. 
lation oTor tho wing longth, Tho scale in these diagrams 
is choson so that tho maxinnim circulation, occurring at tho 
wing roots, is tho same in all casos. She distr i'butions 
for oxtremo high-wing comhinations with iainimuni> induced 
drag aro takom from tho preceding four figures, while those 
for h£gh-wing comhinatlons with conrtanfe circulation, have 
"boon plotted from tho results of the analysis in tho appen- 
dix, .Tho other four cases for mldwing comlainations are^ 

takon. f rom tho results of Lennertz. In those figures, the 
lift distrihu-tions over the fuselage width for the high- 
wing combinations are given by the curves lying below the 
horizontal axis and iaside the vertical lines marking the 
fuselage width. 

As regards the load on the wings, the principal dif- 
. foronce bctwoon tho oxtrom© high-wing (or low-wing) combi- 
nation, and tho midwing combination of tho same span length, 
2b, is that tho formor possess larger wing longth, 
2(b - ain. p tt) . Thus, for a given total load on tho combi- 
nation, and a given span longth, tho load on tho wings is 
higher for tho oxtromo high-wing (or low-wing) combination, 
tho moro so because for those combinations, the lift on tho 
fuBolago is negative so that tho load on tho wings is actu- 
ally greater than the total lift on tho combination. 



Tho fundamental difference between tho extreme high- 
wing (or low-wing) and midwing combinations is that in the 
high- or low-wing eases, the fuselage lift is negative, 
while in the midwing, it is positive. a?his interesting 
result is indicated by the curves shown In figure 19. 
Those curves show the dopendonco of tho ratio of fuselage 
lift to total lift on the combination, that is, Ly/L, on 
tho sonispan b for various types of combination, Tho 
curves, A and B, for midwing combinations are taken from 
the resul^ts of Lennertz; the curves, 0 and D, for ex- 
treme high-wing combinations, are found, respectively, from 
equation (A-15) of the appondix and from the analysis of 
this SQction; tho curve for a combination of internodiato 
wing height with constant circulation has been found from, 
equationi (A-IS) of tho appondix. (Tho portlnont values 
for tho curve 3> aro p tt = 20,5°, sin p tt = 0.350, 
cos p TT = 0,937,) Theso curves tend toward tho value, 
zero, as b increases indefinitely; thoso for tho midwing 
combinations change more slowly than the others. The rea- 
son for this behavior is to be found in equation. (A-16) of 
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the appendix, ffor the midvlng cost'bination, Lp approaches 
a finite linit as becomes infinite, while in the case 

of the extreme high-wing comhination* Lp approaches zero« 

Iho curves for intetbaHato wing heights will cross tho 
axis In general ("but see the following discussion); those 
for the extreme high~wing and midwing comhinaliions do not. 
Also, for the extreme high-wing and midwing comhinations, 
Lp/li is numerically larger in tho case of minimum induced 

drag than in tho case of constant circulation. 

Sho analysis prosonted in tho appendix indicates some 
interesting conclusions in tho case of intornodiato wing 
heights, Tor 1 > cos p tt > 0, tho lift distribution over 

the fuselage is positive over that portions of tho fusolago 
Ixing between tho vortical sections through tho wing roots 
and is negative over tho rost of the fusolago, Tho tran- 
sition in tho loading botwcon these portions occurs "by 
moans of a jump in tho distribution, as shown later in. 
fig:are 23. Although tho distributions for these intermedi- 
ate cases have boon found analytically only for tho case 
of constant circulation (eq.uatiQn8 (A-8) and (A-9) of tho 
appendix), a general consideration of the potential distri« 
bution over the fuselage boundary shows that an exactly 
similar result is obtained for intermediate wing heights 
in tho case of minimum induced drag, lor tho ideal combi- 
nations considorod herein with any listribution of circu- 
lation, there is generally a jump in the loading distribu- 
tion in tho vortical piano through the wing root whoso 
magnitude is P Y r^, where Fj^ is tho circulation at tho 

root, Whothor tho fusolago lift is positive, negative, or 
zero depends on tho relative magnitudes of tho aroas under- 
neath tho separate portions of tho distribution: curves; 
and tho areas, in turn, depend on the wing span, tho wing 
hoight, and tho distribution of circulation alon>g tho wings, 
which may bo constant or be that corresponding to minimum 
induced drag, etc. 

The circumstances under wuich the fuselage lift van- 
ishes are easily determined for the case of constant cir- 
culation. Trom equation (A-12) of tho appendix, let 

Lj. = 2 P V r (sin p TT - rs Vt—) = ° 

^ b + cos P TT'^ 

If this equation is solved for b. 
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■fai = Bin 6 IT, "ba = '■ sin 6 tt (82) 

sin p TT 

The first solution is trivial, for in this case the -wing 
length, 2 ("b - sin p tt) , vanishes so that the total lift 
on the comlJination vanishes as well* Xhe second solu- 
tion is the desired relation. Figure 20 shoTrs tho, dopond- 
once of this "critical somispan" on sin p tt., Tho dashod 
lino of figuro 20 roprosonts tho equation, 

h = sin p TT 

Tho point of intersect iom of the two curves is found from 
tho oq.ua t ion. 

Bin p TT = - ' sin p TT 

sin p TT 

and lies at 

sin p TT = y^/S, or P TT = 45° 

Tho graph shows that for P rr > 45°, "bg < sin; p tt, i,. o,, 

in this range, any wing span necessarily exceeds the crit- 
ical value an^ the fuaelagd lift is necessarily positive,, 
ffor p TT < 45 , the fuselage lift is positive, zero, or 
negative, accordingly as "b > hg , "b = "bg , sin p ti? < "b < "b 

Por tho oxtromo high-wing (or low-wing) combination, 
sin p TT = 0 and "bs is infinite, so that tho fusolago 
lift is negative. In short, for wing heights such that 
P TT > 45°, if tho i7ing span is incroasod, tho fusolago 
lif t ..docroases nunorically hut remains positive; for wing 
heights such that p tt < 45°, tho fuselage lift is nega- 
tive for sufficiently small spas and, as tho span incroasos 
passes through tho value zero whon h = "bg , and thon "be- 
comes positive, Tho second case is illustrated hy tho 
curve E of figuro 19, 

Finally, it should "bo remarked that, if tho comhina- 
tion is regarded as a single structure, tho loading dis- 
trihution ovor tho width of this structure is continuous 
through tho root soctioai, Tho reason for this result is 
that, at tho roots, tho load por unit length of tho wing 
passes from. P V Pj^ to zoro; these discontinuities just 

■balance those in tho fuselage loading. In any case, thoro 
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irill Ijo discontlnuitioe in tlio loading ovor th.o iridtli of 
any particular croBS Boctlom of tlio fusolago at the verti- 
cal sootions that pass through tho roots. 

She analysis prosentod in this eectlon will ho pro- 
foundly modified for coahinations with different fuselage 
cross sections* I'or example, if tho soctions tiro rectangu- 
lar , it can ho readily soon that for tho oxtromo high-xiring 
or lou-vring comhinations, in which the lifting lino is 
tangent to tho upper or the loiror surfacot respoctivoly, 
tho fuselage lift is positive and that portion of tho lift- 
ing lino tangent to tho fuselage ceases to act as a trlng. 
Ihus, tho loading dlstrihutions of those comhinations will 
"bo totally different from those of tho corresponding cem- 
"binationa 'sith circular fuselages, Tho disconnected com- 
"binations with rectangular fuselages will* howover, ho very 
similar to thosQ with circular fuselages. Sxporimonts 
(roforoncoB 9 and 15) also rovoal characteristic diffor,- 
oncos arising from tho shape of tho fuselage cross soction. 
Shoso difforonoos in actual comhinatiens arise from: causes 
au&to different from those doscrihed horca 
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APPENDIX 

EI&H- AKD LOW- WING COHSINAXIONS 
WIIEH OONSIAMT OIHOULAIION 

Again the com'blnatlon 8h.own In figure 2 is considered 
"but now constant circulation oTor the lifting lines is as- 
sumed, Trom the Prandtl theory, the vortex sheet in this 
case dogonoratos into a "•horsoshoo vortex, so that tho 
flow in tho piano at z s 4- <» is that arising from tho 
two vortox filaments trailing in straight linos from the 
wing tips, with the fuselage cross soction as a houndary, 
Eho flow in this plane is ohtainod "by reflecting thoso vor 
tices in the circle, so that the resulting vortex system 
has tho form shown in figure 21 • Tho vortices outside tho 
circle havo tho coordinatos, 

X = ± h, y = COS p TT (A-l) 
and thoso insido havo tho coordinatos. 



c o 



whore 

= TS^ + COS^ P "FT - (A-3) 

In ordor to find tho lift distri'bution over tho fusb- 
lago section, tho relation. 



* 

is used. Tho potential distri'bution over tho fusolago is 
found from olomontary potential theory,, hut caro must "bo 
taken in oxprossing it mathomatically hocause of its multi 
pie values, . The multiple values are avoided here "by the 
introduction of a cut in the z-plane between tho two vor- 
tices outside the circle, . It can bo shown from simple 
considerations that this cut must havo tho form of tho pro 
Joctlon of tho lifting line . on this plane, 7or the case 
considered here, the cut is a horizontal . straight liiie, as 
shown: in figure 21, . 
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Th.e potential field of the vortices outside the cir- 
cle Is: 

= ^ [tan- V -^cos^p n . ,,,-i ^ V: V " ] 



2iT X + (y - COS p tt) - "b 



Vrltlng 



T, »s _ 4.„„ 2b (y COS ^T^) 

X* + (y - cos Ptt) - ^ 



and restricting the tan""^ function to Its principal val- 
ues, -tt/2 < tan"^ Q < Tr/2, the single-valued expressiom 
for this potential is: 

= ~ (x.y), + (y - cos Ptt)® > 
*i = ^ ["■*'^'- ^^'^^J* 2:^+(y-°0B prr) 5 r y> ©os prr > (A-4) 
$1 = ^ [""fT+j^i ('^» x^+(y-cos Ptt)^ 5 h^ I y5 cos §tt 

Shus, for izl this potential is continuous in pass- 

ing through y s cos Ptt> while for )zl < h, it Increas- 
es hy r In passing through this value. 

Similarly writing 

2 ^ (y - oos^N 

■Bi f \ ^ cN. a ^ 
F (x,y) = tan g g- 

s / cos Btt'N "b 



X + 



(j - ^) - 



With the same restriction on the tan"^ function, the po- 
tential of the image vortices can he written as 
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CO s grr 



(A 



r r _, - n1 a./ COS Btt\ ^ "b < cos Btt 

^2 =-2;7[-Tr+raU.y)J.== +(y < ^'y^ — 

Thus, for Ixl > "b/c^, this potential is continuous in 

passing through" y = £SJ^TL, while for Ix' < "b/c^, it 

c 

decreases "by r in passing through this value. In par- 
ticular, this potential is continuous on the circle. 

She total potential is 

$ = $1 + $a 



and its distribution over the fuselage circle, 

+ = 1 

may be written as 

$ = X [tB.n"'- 2b (y - cos Ptt) 

^ 1 + cos Ptt - b - 2y cos prr 

8„ . ££^) 



- tan""^ 



a (cos-* pxT-b^) 

c + % 2y cos Ptt 

c 



where, if the tan~ functions are restricted to their 
principal values, the conditions given in equations (A-4) 
and (A-5) must be employed, This distribution is discon- 
tinuous at the coordinates, x = i sin Tit y = cos tt» 
which correspond to the wing roots, the potential increas. 
ing by the value V , in passing upward through these 
points. Therefore, the potential distribution can bo 
written as: 
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where 



* = ~ [^TT + G(y) j, y > cos Ptt (A-6) 
4 = - TT + G(y)j, y < COS pTT (A-7) 



&(y) = 2 tan 



2y - cos Pn 



In this expression* however, values of the tan""^ function, 
are to he chosen so as to malce G-(y) continuous over the 
circle* Trom the equations, (A-6) and (A-7), the lift 
distrlhution over the fuselage width is found to he: 

PVrf"- ^ — (C^ - 1) V 1 - X** 

dx 



r, ^ * -1 4b (c^ .. 1) ^/ 1 - X^ 1 

I 1 — — tan . ■■ .. .. . . ^. ^ .... .. ^ ^ , 

L " c +1-2 (h -COB Btt) - 4c (1-x ) -1 



Ixl • < Bin p TT (A-8) 

dL™ 1 4h (o® - 1) yi - X® 

w c +l-2Cb -COB p-n) - 4c (1-x^) 

Jx! > sin p TT (A-9) 

The error made hy Lennertz (reference 5) in treating 
this prohlem was apparently caused hy hils failure to sep- 
arate the potential distribution over the fuselage surface 
into distinct parts, so that he obtained equation (A-8) as 
the lift distribution over the entire width, 

These results for the lift distribution are readily 
put into graphical form. Let 



X _x 4b (c^ - 1) y 1 - 



!• (x) = - tan 4 2 2 s a 

^ c +1-2 (b -cos prr) - 4c (1-x ) 

This' function has the form shown in figure 22. Hence, if 
CdLji/dx/pvr) is plotted against x, the curve shown in 

figure 23 is obtained. In the special case of the midwing 
combination (cos P tt = O), this distribution reduces to 
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tHe form obtained Uy Lennertz. shown in figure 2f' J« 
case of the ejctreme high-wing (or low-wing) ??^^^;t*^°Lwn- 
fcos 6 TT = 1). the distrihution reduces to the form shown: 
in figure 25. These last two dlstrihutions appear in fig- 
ures 16 and 18 • 

In order to find the lift on the fuselage, it is un- 
nocossary to integrate tho lift distrihutioa "^JjJ"^^*^ 
has done. Instead, tho total lift on tho comhination is 
found first. Irom oq.uation (26), 



I, = P V H 



Writing 



.P./ 

00 



f(z) dz 



Z± = "b + i cos P TT 



the flow function in this case is: 



f(z) = 



r 

2iTi 



Z — Zt _ 

log =r- + log 

z + Zj. 



s — 



_ r 



2TTi 



Z, + Z. 



log (l - ^^441) + log X + — s-^ 



z - 

c 



Expand in descending powers of (z + "zi) and (z - 



f'(z) = 



2rTi 



2Tri 



z, + z, (z,+Zj)/c^ 

_ — i i. + = ■ ■ ■ + • • • 

Z + Zl 



2h 



-1 , ■ 1/c . 

^. — - t 1111 ■!■ + • • * 



Z + Zi 



Zi 



- ^ 



as 



Zi + Zi =! 213 
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Tlien 

- P T r 21. (1-^) 

■Shis expression for the lift is based on nondimeuBlonal 
lengths; the dimensional relation le 

L=PVrE21) (1- ^) (A-10) 

vrhero E is the fuselage radius. Thus, for the same cir- 
culation and span length, the extreme high-irlng (or low- 
wing) comhination has greater lift than the mldvlng comhi- 
natiost although for practical values of h, the differ- 
ence is negligible; and, in either case, the lift does not 
differ greatly from that of an Isolated wing of the same 
span, 

She lift on the wings is: 

Lir = P V r H 2 (h - sin p tt) (A-11) 

Hence, 

Li. = L-Lff = 2 PVTR ^sin p tt - (A-12) 



and 



- a 1> 
sm p IT - 72 sf 



1, g ^ 



^. + P (A- 13) 



"b + cos p TT 

for vx.e mldwing combination, 

sin p TT = 1, cos p TT = O 



and 



^7 



L b + 1 



(A-14) 



Ihis result is shown graphically by the curve B of figure 
19« ffor the extreme high-wing (or low-winji) combination. 



sin p TT = 0, cos p ir = ± 1 
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so that 



1 



(A-15) 



Shis result is shown "by the carve 



of figure 19," 



She limiting case of infinite wing span is of special 
interest. For this case, equation (A-12) reduces to 



Lj = 2 P T r H sin p tt 



(A-16) 



Thus, Dp vanishes for the extreme high-wing and low-T7ing 
combinations with infinite wing span, while for the midwlng 
combination. 



Lj. = 2 p V r R 



(A- 17) 



The general result, equation (A-16), signifies that for in- 
finite span tho lift on the combination is just the same 
as if the wings were continuous through the fuselage and 
the fuselage removed. The particular result, equation 
(A-17), was derived by Lonnertz (reforonco 5) by an in- 
volved mathematical analysis and has since been verified 
by experiment (reference 11). 

The analysis presented here applies as well to dis- 
connected combinations, i.e., those in which the wings do 
not intersect the fusolago but lie at some distance above 
or below it« If tho nondimonsional height of tho lifting 
lino from the fusolago axis is callod h, it is roadily 
soon from equations (A-10) and (A-13) that, in this case. 



and 



so that 



L = 2 P 7 r R b (l - 

Lj = 2P7rRb(-.^^-L_) 



L ^ 



+ h - 1 



(A-18) 
(A-19) 

(A- 30) 



Thus, for such combinations, is always negative 

and vanishes for infinito wing span* 

This thoorotical result agrees poorly with tho measure 
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NACA a>echnlcal Hote Ho. 812 



41 



ments of the foroes acting on the separate memherB of dis- 
connected comhinations. In the tests of reference 12 it 
was found that the Interference lift force on the fuselages 
of disconnooted high-wing comhi^ations wore predominantly 
positiTs, She reason for this poor agreement appears to 
lie in the oreationi of a low-pressure region between wing 
and fuselage hy means of a venturi effect arising from the 
finite profile of the airfoil and the curvature of the 
fusolage in side elevation. These quantities do not appear 
in the first approximation of the theory used here. She 
experiments also reveal (reference 12) tho presence of an 
Interforenco lift forco acting on tho wlng« of which no in- 
dication! appears in tho theory. 

Finally, the lifting-line theory yields another inter- 
esting result for comhinations with constant circulation. 
It follows directly from the last analysis that the lift 
distribution over the entire width of the combination, 
that is, over both fuselage and wings, depends only on the 
po-sitions of tho wing roots and wing tips and is ontiroly 
Independent of the front elevation of the wings. 

Other results of the application of the lifting-line 
theory to combinations with finite fuselages have been 
obtained by 7. Tandrey (reference 16). 
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Figs. 6,7,11,12 
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Y Figs. 21,22,23,24,25 
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